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We study the classical dynamics of the Abelian Higgs model employing an asymptotic multiscale 
expansion method, which uses the ratio of the Higgs to the gauge field amplitudes as a small 
parameter. We derive an effective nonlinear Schrodinger equation for the gauge field, and a linear 
equation for the scalar field containing the gauge field as a nonlinear source. This equation is 
used to predict the existence of oscillons and oscillating kinks for certain regimes of the ratio of 
the Higgs to the gauge field masses. Results of numerical simulations are found to be in very 
good agreement with the analytical findings, and show that the oscillons are robust, while kinks 
are unstable. It is also demonstrated that oscillons emerge spontaneously as a result of the onset 
of the modulational instability of plane wave solutions of the model. Connections of the obtained 
solutions with the phenomenology of superconductors is discussed. 
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1. Introduction 

Soliton solutions of field theoritical models play a significant role in the description of physical 
phenonena occuring in a wide class of systems ranging from high energy physics to cold atoms and 
granular media [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. 
In particular classical solutions of the mean field Gor’kov - Eliashberg - Landau - Ginzburg theory 
are important for the phenomenology of superconductivity [26, 27, 28, 29] where the emerging 
effective field theory is the Abelian-Higgs model. Several solutions of the latter are found [30, 31, 
32]. 

In the present work [33], we analyze the classical dynamics of the Abelian-Higgs model in 
(1 -|- l)-dimensions by means of a multiscale expansion method [35]. We restrict our analysis on 
the case of small fluctuations of the condensate about its vacuum expectation value (vev). This 
is realized by choosing the amplitude of the Higgs field to be an order of magnitude smaller than 
that of the gauge field. In this limit the dynamics simplify considerably and it is found that the 
scalar field performs asymmetric oscillations around the classical vacuum. Such a scenario occurs 
naturally considering the model just after the symmetry breaking i.e. close to the critical point. 
Then the minimum of the scalar field potential is very flat and the asymmetric cubic term is strong 
leading to an asymmetric shape of the potential around it. However, results are also found in 
the case where the gauge and the scalar field amplitudes are of the same order [34]. This scenario 
corresponds to a strong breaking of the underlying gauge symmetry, which is far beyond the related 
critical point. In this case the minimum of the potential occurs at the bottom of a deep well, while 
the potential shape is almost symmetric around it since the quadratic term dominates. 

Employing the method of multiple scales (see also Ref. [36]), the original nonlinear coupled 
field equations are reduced to an effective nonlinear Schrddinger (NLS) equation for the gauge field, 
and a linear equation for the Higgs field containing the gauge field as a source. These equations are 
analytically solved giving two types of localized solutions for the gauge field, namely in the form 
of oscillons and oscillating kinks. Subsequently, we numerically integrate the original equations 
of motion using, as initial conditions, the analytically found solutions. We find that the analytical 
predictions are in excellent agreement with the numerical findings for a large range of values of 
the parameters involved. We then discuss connections between the Abelian-Higgs model solutions 
and the phenomenology of superconducting materials, describing the temporal fluctuations of the 
condensate in a one-dimensional (ID) Josephson junction [37] due to the presence of oscillating 
(in time) magnetic and electric fields. 

2. The model and its analytical consideration 

2.1 Formulation and Equations of Motion 

The t7(l)-Higgs field dynamics is described by the Lagrangian: 

^ = + (D^(j))*(D^^(j))-E(<I>*d>), (2.1) 

where F^y is the 17(1) field strength tensor, + ieA^ is the covariant derivative, and 

E(<I)*d)) = -|- A(<I>*d>)^ (A > 0) is the Higgs self-interaction potential (asterisk denotes 
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complex conjugate). In the broken phase, < 0, a vev V js/l of the Higgs field arises classically: 

= —pp-jX. We will focus on the dynamics of this system assuming that the Higgs field fluclu- 
afes slighfly around ifs vev. In fhis case, we expand fhe field <I> as: <I> = ^(t) +^) and obfain fhe 
following equations of motion for and H'. 


(n + ni^A^ — d^idyA^^-\-'2.e^vHA^-\- —0, 

(□ + m]j)H + 3XvH^+ Xh'^- e^A^A^^{v+H) = Q, 


( 2 . 2 ) 

(2.3) 


where m\ = and = 2Au^. We simplify Eqs. (2.2)-(2.3) choosing fhe field represenfafion 
Aq = Ai = A 3 = 0, and A 2 = A / 0. Due fo fhe facf fhaf we are interesfed in ID sellings, fhe non 
vanishing A 2 field is a function solely of x and t. As a consequence, fhe Lorenfz condition dyA'' = 0 
is frivially fulfilled and we are leff wifh a coupled sysfem of equations for fhe gauge field A(x,t) and 
fhe Higgs field. Furlhermore, we wrife Eqs. (2.2)-(2.3) in a dimensionless form by rescaling fhe 
fields as: A —)■ (in/^/e)A, H — {mA/e)H, and space-lime coordinafes as: jc —)■ x/iua and t — t/niA- 
Thus we end up wifh fhe following equafions of molion: 

(□ + 1 )A + 2AH + H^A = 0, (2.4) 

(□ + q'^)H+^q'^H^ + ^q'^H^+HA^+A^ = 0, (2.5) 

where parameter q = mu/mA is assumed lo be of order ^(1). The Hamiltonian corresponding to 
the above equations of motion is: 

i [{d,Af + {dAf + {dtHf + [d^Hf] + E, (2.6) 

where indices denote partial derivatives with respect to x and t, and the potential V is given by: 

V = —h^ + —h'^ + —h^ + -h'^a^ + ha^ + -a^. 

8 2 2 2 2 

(2.7) 


Notice that V exhibits a single minimum at {A,H) = (0,0). 


2.2 Multiscale expansion and the NLS equation 

Considering a small fluctuating field H, we may employ a perlurbalion scheme, which uses a 
formal small paramefer 0 < £ <C 1 defined by fhe rafio of fhe amplilude of fhe Higgs field fo fhe 
amplilude of fhe gauge field, i.e., e ~ H/A. In particular, we will employ a multiscale expansion 
melhod [35], assuming lhal fhe Higgs and gauge fields depend on fhe sef of independenf variables 
xq = X, xi = ex, X 2 = s^x,... and to = t,t\= et, t 2 = e^t ,.... Accordingly, fhe partial derivative 

operators are given (via the chain rule) by dx = 5^0 + + ■ ■ ■, dt = dt^ + £5?, +_Furthermore, 

taking into regard that the fields should be expanded around fhe frivial solufion of Eqs. (2.4)-(2.5), 
as well as ~ £A as per our assumpfion, we infroduce fhe following asympfolic expansions for A 
and H: 

A = £A^'^ + e^A^^^ +..., 

H = + + ..., 
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where and (7 = 0,1,...) denote the fields at the order Substituting Eqs. (2.8)-(2.9) 

into Eqs. (2.4)-(2.5), and using the variables vq, to> • •> we obtain the following system of 


equations up to the third-order in the small parameter epsilon 

^(£) : (□o + 1)a(') =0, (2.10) 

^(£2): (□o + i)aP)+2(5,„6),, =0, (2.11) 

(□o + ^2)//(^^+a(i)2 = o, (2.12) 

^(£3): (□o + 1)a(3)+2(5,„4 

+ (ni+25,„4-24„4+2//P))A(') =0. (2.13) 

The solution of Eq. (2.10) at the first order is: 

A^'^ = m(vi,V 2 , . .. ,t2,t3 ■ ■ +C.C., (2.14) 


where u is an unknown function, “c.c.” stands for complex conjugate, while 6 = Kx — Q.t, and 
the frequency and wavenumber K are connected through the dispersion relation Q? = + 1. 

The solvability condition for the second-order equation (2.11), is (dt^dti — = 0 ^ which 

satisfied in fhe frame moving wifh fhe group velocify Vg = dQ./dK, so fhaf u depends only on fhe 
variable xi =xi — X)gt\ and t 2 - Then, fhe field A^^) is taken to be of the same functional form as A^*) 
[cf. Eq. (2.14)]. Taking into account the above results, Eq. (2.12) for the Higgs field becomes: 

=B{b\u\^ + u^e-^^+ c.c), (2.15) 


where B = 1/(4 — q^) and b = —IjBq^-, here, if is assumed fhaf q ^0 and ^ / 2, as for fhese 
limifing values fhe field //( 2 ) (which is assumed fo be of i^(l) for fhe perfurbafion expansion fo 
be valid) becomes infinilely large. Einally, af fhe order ^(£^), we firsl nofe fhaf fhe second ferm 
vanishes in fhe frame moving wifh u. Thus, fhe solvabilify condifion (obfained —as before— by 
requiring fhaf fhe secular pari is zero) becomes: (Di + ldtQdt^ — Idx^dx^ + = 0. In fhe 

following, for simplicify, we will consider fhe zero momenfum case (K = 0); hence, Vg = 0, and 
also jci =xi, i.e., fhe field A 1 is af resl. In fhis case, we derive fhe following NES equafion for fhe 
unknown funcfion u{x\,t 2 )- 



where fhe parameter s is given by 
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2 1 

q^~^ q^ —4 


(2.16) 


(2.17) 


The above NES equafion possesses exacl localized solufions of fwo differenf fypes, depending 
on fhe sign of s: for 5 > 0, fhe solutions are sech-shaped solifons, while for 5 < 0 fhe solutions are 
fanh-shaped kinks. These solutions have been extensively sfudied in a variety of physical confexfs 
including nonlinear optics [38], wafer waves [39], atomic Bose-Einstein condensates [25], and so 

*The second term is “secular”, i.e., is in resonance with the first term which becomes == 6exp(i6), thus leading to 
blow-up of the solution for t +<». 
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on. Here, the type of solution, i.e., the sign of the nonlinear term in Eq. (2.16), depends solely on 
the parameter q. 

In particular, for ^ < 1.63 and q > 2, we obtain 5 > 0 [attractive (focusing) nonlinearity] and 
the soliton solution of Eq. (2.16) has the form: 

u = uosech[^/suoXl) , (2.18) 


where uq is a free parameter [considered to be of order ^(1)] characterizing the amplitude of the 
soliton, while the soliton frequency is ft) = —{\/2)u^s. Accordingly, the approximate solutions of 
Eqs. (2.4)-(2.5) for A and H can be expressed in terms of the original variables x and t as follows: 


A « 2£ftosech [euqs/sx) cos 
H K, (£Mo)^®sech^ (£MoV^-^) \ b + 2 


1 - ]^{euQfs ) t 


cos 


2 (1 --{euofs ) t 


(2.19) 

( 2 . 20 ) 


characterized by the single free parameter euq. The above approximate solutions for the gauge 
field A and the Higgs field H are localized in space (decaying for |x| —)• °°) and are oscillating in 
time; thus, they correspond to oscillons (alias breathers). At the leading-order in £, the oscillation 
frequencies for A and H are given by: ftjq ~ 1 and (Oh ~ 2 (in units of ma) respectively. 

On the other hand, for 1.63 < q < 2, we obtain s <0 [repulsive (defocusing) nonlinearity], so 
that Eq. (2.16) possesses a kink solution of the following form: 

ft = ftotanh (2.21) 

where fto is a ^(1) free parameter and (O = Uq\s\. In this case, the approximate solutions for the 
gauge and Higgs fields read: 

A « 2£ftofanh (^s/\i\Euox^ cos [(l + (£fto)^|5|) t] , (2.22) 

H ^ (£fto)^Bfanh^ ^Y^£ftox^ {^ + 2cos [2 (l+ £^ftQ|5'|) t] } . (2.23) 

These solutions correspond to oscillating kinks, for both fields, which are zero exactly at the core of 
the kink and acquire non vanishing values at \x\ —)• oo. Again the fields oscillate with the frequencies 
ftjq ~ 1 and (Oh ~ 2. 


3. Numerical results 

To further elaborate on our analytical findings, in this Section we will present results stemming 
from numerical integration of Eqs. (2.4)-(2.5). Our aim is to check the validity of our analytical 
predictions, namely the existence of oscillons and oscillating kinks, as well as study numerically 
their stability. The equations of motion are integrated by using a fourth-order Runge-Kutta time 
integrator and a pseudo-spectral method for the space derivatives [40]. The lattice spacing used 
was fixed to Ax = 0.2, the time step At = 10^^, and the total length of the lattice L = 400. In all 
simulations we used, as an initial condition, the approximate solutions of Eqs. (2.19)-(2.20) for the 
oscillons and of Eqs. (2.22)-(2.23) for the oscillating kinks for ? = 0; in all cases we have fixed fhe 
free parameter value to fto = 1, and the ratio of the two fields amplitudes to £ = 0.1. 
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Figure 1: (Color online) Top panel: Evolution of the gauge field A as a function of x and t. Bottom panel: 
Profile snapshots of the field at different instants of its breathing motion. Parameter values used are £ = 0.1, 
^=1.5 and the oscillation period is T 4 = 2n. 




X 


Figure 2: (Color online) Same as in Fig. 1, but for the Higgs held H. Notice that the oscilation period of the 
Higgs held is Th = K. 


3.1 Oscillons 

In the top panels of Fig. 1 and Fig. 2 a contour plot showing the evolution of an oscillon 
solution is shown (for ^ = 1.5) at the end of the integration time t = 5 x 10'*. Note that the periods 
of the two fields are TA = 2n/(OA = 2n and Te = Ih/coh = tt, and thus 3 oscillations for the gauge 
field and 6 oscillafions for fhe H field are depicfed. In fhe bottom panels of Figs. 1 and 2, we show 
fhe profiles of A and H , during fhe interval of one period of exhibiting fheir breafhing mofion. 
While fhe gauge field A performs symmefric oscillations, fhe Higgs field oscillates asymmefrically 
wifh respecf fo ifs equilibrium poinf. This is due fo fhe consfanf b in fhe solufion of Eq. (2.20). The 
evolution of an oscillon for q = 2.5 (lying in fhe second region where oscillons exisf) is also shown 
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Figure 3: (Color online) Same as Fig. 1, for parameter values e = 0.1 and q = 2.5. 

in the top panels of Fig. 3 and Fig. 4. The breathers become more localized in this region, as can 
be seen from their snapshots (bottom panels of Figs. 3-4). 

In order to demonstrate the localization of the energy of the oscillons, in the top and middle 
panels of Fig. 5, we show a contour of the total Hamiltonian density [cf. Eq. (2.6)] for ^ = 1.5. The 
top panel, corresponds to an initial condition using only Eqs. (2.19)-(2.20), while in the middle 
panel we have also added a random Gaussian noise to the initial condition, as large as 10% of 
the oscillon’s amplitude (the rest of the parameters are as in Eig. 1). We observe that the energy 
density remains localized during the numerical integration in both cases (with and without noise). 
We would like to stress here that, the latter result, shows that oscillons are robust under the effect 
of a random noise, and this was also confirmed for different values of q in the domain of existence 
of the oscillon. Eor completeness, the total energy E{t) = normalized to its initial 

value £"(0), is also shown (dashed dotted line)in the bottom panel of Eig. 5. In fact the energy 
fluctuations, defined as AE = E — E{0), as shown in Eig. 5 [for ^ = 1.5 (top solid line) and q = 2.5 
(bottom dashed line)] are of the order of lO^'*. 

3.2 Spontaneous oscillon formation 

In this section we show that oscillons can emerge spontaneously, through the mechanism of 
modulation instability (for more details see Ref. [41]). The latter concerns the instability of plane 
wave solutions of the NES Eq. (2.16), of the form 

u{xi,t2) = UQe,\p[i{kx\ — cug)], (3.1) 

under small perturbations. To briefly describe the emergence of this instability, we consider the 
following ansatz: 

u = {uq + VF) exp [/ {kx\ — (Ot 2 ) + (0)], (3.2) 

where the amplitude and phase perturbations W and 0 are given by: 

W = £Wbexp[/(2xi — Pt2)] +C.C., (3.3) 
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Figure 4: (Color online) Same as Fig. 2, with parameter values e —Q.l,q — 2.5. 



Figure 5: (Color online) Top and middle panel: contour plot showing the energy density J^(t), as a function 
of time, for ^ = 1.5. For the middle panel the initial condition was perturbed by a small random noise with 
an amplitude 10% of the oscillon’s amplitude. Bottom panel: the total energy normalized to its initial value 
at f = 0 dashed dotted line [upper (black)] . The energy difference AE, for ^ = 1.5 solid [middle (blue)] line 
and q = 2.5 dashed [bottom (red)] line. 


0 = £0oexp [/(2xi — Pf2)] +C.C., (3.4) 

with Wo and 0o being constants, and e being a formal small parameter. Substituting Eqs. (3.2)- 
(3.4) into Eq. (2.16) it is found that, for k = 0, the frequency P and the wavenumber Q of the 
perturbations, obey the dispersion relation: 

P^ = {l/2)\Q\WQ\^/2-2s\m\^). (3.5) 

It is evident from Eq. (3.5) that for 5 > 0 there exists an instability band for wavenumbers < Au^s, 
where P becomes complex. Note that this instability can only emerge in the region where oscillons 
exist. When this instability manifests itself, the exponential growth of the perturbations leads to the 
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Figure 6: (Color online) Top panel; Contourplot showing the energy functional J^{t), for q= 1.5, 2 = 0.1, 
Mo = 1, e = 0.1 for an initial condition corresponding to an unstable plane wave. Bottom panel: Prohle of 
the gauge held A (.r), (solid [black] line), at f = 13040 after the instability has set in. Dashed (blue and red) 
lines, correspond to the the solutions of Eq. (2.19), htted with the parameter mq at the peak of the oscillons 
(indicated by the dotted [grey] lines at ;<: = 75 and x = 150). 


generation of localized excitations, which are identified as the oscillons described in Eqs. (2.19)- 

( 2 . 20 ). 

To illustrate the above, we have numerically integrated Eqs. (2.4)-(2.5), with an initial condi¬ 
tion corresponding to the plane wave of Eq. (3.1), perturbed as in Eq. (3.2), with e = 0.1, Wq = 1, 
2 = 0.1 (inside the instability band), and 0o = 0; the rest of the parameters used are ^ = 1.5, 
e = 0.1, Mo = 1 and k = 0. In the top panel of Eig. 6, we show a contour plot of the energy den¬ 
sity, It is observed that, at t ~ 10,000, localization of energy is observed due to the onset 

of the modulation instability. This localization is due to the fact that harmonics of the unstable 
wavenumber Q of the perturbation are generated, which deform the plane wave and lead to the for¬ 
mation of localized entities. In fact, the latter are eventually reformed into oscillons, as is clearly 
observed in the bottom panel of Eig. 6: there, we show the profile of fhe gauge field A{x) [solid 
(black) line] af t = 13040, and we idenlify af leasf fwo well formed oscillons located af x = 75 and 
X = 150. The laffer are found fo be in a very good agreemenf wifh fhe analytical profile depicfed 
by fhe dashed (red and blue) lines. To perform fhe tiffing, we plotted fhe approximate solufions of 
Eq. (2.19), using fhe value af fhe peak of fhe oscillons, in order fo identify fhe oscillon amplifude 
Mq. Then, fhe solufion was also displaced by a consfanf faclor, in order fo mafch fhe cenfer of fhe 
oscillon. Addifionally, time was sef fo zero (t = 0) and fhe fitting was made af fhe beginning of fhe 
period of oscillafion for fhe particular oscillons. The very good agreemenf befween fhe analyfical 
and numerical field profiles highlighfs fhe facf fhaf fhe oscillons considered in fhis secfion can be 
generafed sponfaneously via fhe modulafional insfabilify mechanism. 

3.3 Oscillating kinks 

We now proceed wifh fhe numerical sfudy of solutions corresponding fo fhe kinks of fhe NES 
Eq. (2.16). Such solufions in fhe form of Eqs. (2.22)-(2.23) are expected fo exisf in fhe parameter 


9 











Oscillons and oscillating kinks in the Abelian-Higgs model 


Charilaos Tsagkarakis 



region 1.63 < ^ < 2. Performing the same procedure as in the previous section, we numerically 
integrate the equations of motion and study the relevant dynamics. In the top panels of Figs. 7 
and 8, the contour plot of a kink solution is plotted, showing the evolution of both fields A and 
H, during the interval of three and six periods, respectively (for ^ = 1.8). In addition, the profiles 
of bofh fields are shown in fhe bottom panels of Figs. 7- 8. The oscillafing kinks were found fo 
be unsfable for all values of fhe paramefer q wifhin their region of existence. The instability is 
manifested by an abrupt deformation of the kink solution, even near its core, characterized by the 
inverse width Id, (where U = l-\/\s\euQ). In Fig. 9, the top panel shows the profiles of fhe gauge 
field af fhe beginning {t = 8) [solid (blue) line], and af fhe end of fhe infegrafion {t = 1200) [solid 
(red) line]. While fhe shape of fhe kink is more or less preserved, if is clear fhaf fhe solufion profile 
has been significanfly deformed. More imporfanfly, fhe Higgs field H, shown in fhe boffom panel, 
nol only has been disforfed buf if has also become an order of magnifude larger fhan ifs inifial 
amplifude. Thus, if can be concluded fhaf oscillafing kinks are unsfable. 


4. Discussion and conclusions 

In fhis work we have presenfed a class of classical solufions of the ID Abelian Higgs model 
obtained with the use of a multiscale expansion method. The key assumption in our treatment 
is that the scalar field amplifude is much smaller fhan thaf of the gauge field; fhe relevanf ratio 
was fhen used as a formal small paramefer in the perturbation expansion. We have shown that the 
equations of motion can be reduced to a NFS equation; by means of the latter, localized solutions 
in the form of oscillons and oscillating kinks were derived. Results by numerical integration of the 
original equations of motion where found to be in very good agreement with the analytical findings. 
In addition, we have numerically sfudied sfabilify of fhe solufions againsf a Gaussian noise wifh 
amplifude up fo 10% of fhe Higgs field amplifude, and found fhaf fhe oscillons remain robusf. The 
robusfness of oscillons is also complemenfed by more recenf findings [34] where if is indicafed 
fhaf only oscillon solufions for fhe gauge and fhe Higgs field are long-lived. This leads us fo fhe 
conclusion fhaf oscillons dominafe in fhe solution space of fhe Abelian-Higgs model. 
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Figure 9: (Color online) Top panel: profile snapshots of the gauge field at the beginning [solid (blue) line], 
and at the end of the simulation [solid (red) line], for q= 1.8. Bottom panel: same as in top panel but for the 
Higgs field. 


It is also relevant to discuss the possible connection of the presented solutions to the physics of 
superconductors. One could, in principle, write down the form of the magnetic and electric fields 
originating from the gauge field A which was chosen to be in the y direction: B{x,t) = dxA{x,t)z, 
and S{x,t) = —dtA{x,t)y. 

In the case of the oscillon solutions given in Eqs. (2.19)-(2.20), the above equations describe 
the electric field in the y direction, which produces a magnetic field in the z direction; both fields are 
localized around the origin of the a: axis - cf. Fig. 10, where the profiles of the fields are shown. This 
can be thought of as a configuration, describing Superconductor-Normal metal-Superconductor 
(SNS) Josephson junction [37], where two superconductors are linked by a thin normal conductor 
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Figure 10: (Color online) Top panel; profile snapshots of the electric field S'y{x) for q= 1.5. Bottom panel; 
profile snapshots of the magnetic field B^{x) for the same q. 

placed at the origin. Then, our solutions describe a condensate (j) = {v + H{x,t))/V2 that oscil¬ 
lates around its vev near the origin, and acquires its vev value when entering the superconductors. 
Accordingly, the magnetic field is shown to oscillate inside the normal conductor, but vanishes 
exponentially inside the superconductors as per the Meissner effect. 

Our approach not only reveals a new class of solutions of the Abelian-Higgs model, but also 
dictates a straightforward general strategy for the search of non-trivial dynamics in models in¬ 
volving classical fields with nonlinear interactions. In particular the dynamics of such models is 
governed by the nonlinear Schrddinger equation. It is an interesting perspective to determine the 
impact of these solutions on the thermal and quantum behaviour of the involved fields. Such studies 
is an interesting theme for future work. 

Acknowledgments. Illuminating discussions with L. P. Gork’ov are kindly acknowledged. 
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